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Abstract. 

We analyze the validity of perturbative renormalization group estimates obtained 
within the fixed dimension approach of frustrated magnets. We reconsider the 
resummed five-loop /3-functions obtained within the minimal subtraction scheme 
without e-cxpansion for both frustrated magnets and the well-controlled ferromagnetic 
systems with a cubic anisotropy. Analyzing the convergence properties of the critical 
exponents in these two cases we find that the fixed point supposed to control the second 
order phase transition of frustrated magnets is very likely an unphysical one. This is 
supported by its non-Gaussian character at the upper critical dimension d — 4. Our 
work confirms the weak first order nature of the phase transition occuring at three 
dimensions and provides elements towards a unified picture of all existing theoretical 
approaches to frustrated magnets. 
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1. Introduction. 

Although undoubtedly successful to describe the critical behavior of 0(N)-like 
models, perturbative field theory is still unable to provide a clear, non-controversial 
understanding of the physics of certain more complex models among which are the 
famous Heisenberg or XY frustrated magnets (see [1] and references therein). At the 
core of the problem is that different kinds of perturbative approaches, performed up 
to five- or six-loop order, lead to contradictory results: in dimension d = 3, first order 
phase transitions are predicted within the e (or pseudo-e)- expansion [21 |3l H] whereas 
a second order transition is found in the fixed- dimension (FD) perturbative approaches 
performed either in the minimal-subtraction (MS) scheme without e-expansion [5] or in 
the massive scheme (6j. In fact, FD results for frustrated magnets are neither supported 
by experiments nor by Monte Carlo simulations pfl [7J El E2, HD] (see however [5] where 
a scaling behavior is found). They also disagree with the results obtained from the 
non-perturbative renormalization group (NPRG) approach [U |TT] that predicts (weak) 
first order phase transitions in d = 3 in agreement with the e-expansion analysis. 

In this article we shed light on the discrepancies encountered in perturbative 
approaches to frustrated magnets by showing that the FD approaches lead to dubious 
predictions as for the critical physics in d = 3. Our key-point relies on the very nature 
of the FD perturbative approach and is easy to grasp already for the simplest — O(N) 
- model. In this case, the (non-resummed) renormalization group (RG) /5-function 
at L loops is a polynomial in its coupling constant u of order L + 1. Thus, it admits 
L + 1 roots u* , (3{u*) = 0, that are either real or complex. Within the e-expansion, 
when one solves the fixed point (FP) equation in successive orders in e = 4 — d, the 
only non-trivial FP retained is by definition such that u* ~ e. On the contrary, in the 
FD approaches, when one directly (analytically or numerically) solves the non-linear 
FP equation at fixed d (fixed e) no real root can be a priori discarded. As a result, the 
generic situation is that the number of FPs as well as their stability vary with the order 
L: at a given order, there can exist several real and stable FPs or none instead of a 
single one. This artefact of the FD approach is already known and was first noticed in 
the massive scheme in d = 3 |12j . The way to cope with it is also known: resumming the 
perturbative expansion of (3{u) (see e.g. [13]) is supposed both to restore the non-trivial 
Wilson-Fisher FP and to suppress the non-physical or "spurious" roots. This is indeed 
what occurs for the O(N) model for which the FP analysis performed on the resummed 
/3-function of FD approaches enables to discriminate between physical and "spurious" 
FPs. However this ability of the resummation procedures to eradicate spurious solutions 
of the FD approach has never been questioned and, a fortiori, evaluated in the context 
of more complex models and, in particular, for models with several coupling constants. 
We precisely argue, in this article, that the situation is very different for frustrated 
magnets and probably for several other models. Indeed, considering the /3-functions 
derived at five loops in the MS scheme and using a standard resummation procedure [5] 
we show that the FP found in d = 3 without expanding in e, although it persists after 
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resummation, is in fact spurious. 

Our conclusion is based on several facts: (i) the critical exponents computed at 
the FP supposed to control the second order behavior of frustrated magnets display 
bad convergence properties with the order of computation in the controversial cases of 
XY and Heisenberg spin systems (ii) when analyzed with the same FD approach, the 
field theory relevant to ferromagnetic systems with cubic anisotropy displays a similar 
FP — having no counterpart within the e-expansion — in contradiction with its well 
established critical physics. The critical exponents computed at this supernumerary FP 
display the same bad convergence properties with the loop order as in the case of XY 
and Heisenberg frustrated magnets (iii) the coordinates (u^u^) of the attractive FP 
found in the FD approach of frustrated magnets are multivalued functions of (d, N) - 
N being the number of spin components — because of the existence of a topological 
singularity in the mapping between (d, N) and the FP coordinates (u\(d, N), u^d, N)). 
This singularity provides strong indications of the existence of pathologies in the RG 
equations obtained at fixed dimension (iv) finally, we provide strong arguments showing 
that the supernumerary FPs found in the frustrated and cubic models survive in the 
upper critical dimension d = 4 where they are found to be non-Gaussian, a behavior 
deeply connected with the existence of the above mentioned topological singularity. 
Being given the present state of knowledge of 4 -like theories that are very likely 
trivial in d — 4, this fact confirms the serious doubts on the actual existence of these 
supernumerary FPs. 

2. Resummation method. 

To investigate the five-loop /3 functions derived in the MS scheme we have to resum 
them. Before discussing the case of a series of two coupling constants relevant for 
frustrated magnets or ferromagnetic models with cubic anisotropy we recall, for the 
sake of clarity, the main steps needed to resum a series of one coupling constant u as 
well as the underlying hypotheses (see for a review [II]). Let us consider a series 



and is supposed to converge, in the complex plane, inside a circle of radius 1/a, where 
u = — 1/a is the singularity of B{u) closest to the origin. Then, using this definition as 
well as T [n + b + 1] = J °° t n+b e~ l dt, one can rewrite 




(1) 



n 



where the coefficients supposed to grow as n\. 

The Borel-Leroy sum associated with f(u) is given by: 




(2) 




(3) 
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and, interchanging summation and integration, one can define the Borel transform of / 
as: 

POO 

f B (u) = dt e~* t b B{ut) . (4) 
Jo 

In order to perform the integral in (jlj) on the whole real positive semi-axis one has to 
find an analytic continuation of B{t). Several methods can be used, Pade approximants 
for instance. However, it is generally believed that the use of a conformal mapping is 
more efficient since it makes use of the convergence properties of the Borel sum. Under 
the assumption that all the singularities of B(u) lie on the negative real axis and that 
the Borel-Leroy sum is analytic in the whole complex plane except for the cut extending 
from — 1/a to — oo, one can perform the change of variable: 



VI + au - 1 A uj 
uj{u) = . -<=>- uiuj) = -— — (5) 

that maps the complex w-plane cut from u = —1/a to — oo onto the unit circle in the 
w-plane such that the singularities of B(u) lying on the negative axis now lie on the 
boundary of the circle \w\ = 1. The resulting expression B(u(u>)) has a convergent 
Taylor expansion within the unit circle |u;| < 1 and can be rewritten: 

S(«(w)) = ^^(0,6) [u(u)] n (6) 

n 

where the coefficients d n (a,b) are computed so that the re-expansion of the r.h.s. of 
(EJ) in powers of u coincides with that of (0Q). One obtains through (jHJ) an analytic 
continuation of B(u) in the whole u cut-plane so that a resummed expression of the 
series / can be written: 



/*(«) 



poo 

Vrf n (a,6) / dt e-H h [u{ut)) n . (7) 

In practice it is interesting to generalize the expression (JJJ) by introducing [15] the 
expression 

f R (u) = Y,dn{oc,a,b) f°°dt e-H b (8) 

whose meaning will be explained just below. 

If an infinite number of terms of the series fn(u) were known, expression dSJ) would 
be independent of the parameters a and b and a. However when only a finite number 
of terms are known, Jr(u) acquires a dependence on them. In principle, the parameters 
a and b are fixed by the large order behavior of the series: 

a n ^ oo ~(-a ) n n\n b0 (9) 

which leads to a = a and b ~ b + 3/2 [16] while a is determined by the strong coupling 
behavior of the initial series: 

f(u -> oo) ~ u aa/2 (10) 
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which can be imposed at any order of the expansion by choosing a = olq. The common 
assumption is that the above choice of a, b and a improve the convergence of the 
resummation procedure since it encodes exact results. 

Let us however emphasize that, often, only a is known and that the other 
parameters, a and b, are considered either as free (as for instance in [5]) or variational 
(for instance in [T7] where a is determined by optimizing the apparent convergence of 
the series). In any case, the choice of value of a, a and b must be validated a posteriori 
by checking that a small change of their value does not yield strong variations of the 
quantities under study. Such variations would clearly indicate that one has chosen an 
unstable, parameter-dependent calculation procedure or that one has not computed the 
quantities under study at a sufficiently high order of perturbation theory to consider 
them as converged. In the following, we shall employ this "stability criterion" to validate 
- or invalidate - the results obtained by means of the FD perturbative approach to 
frustrated magnets. 

In the context of frustrated magnets, the above described resummation procedure 
must be extended to several (two) coupling constants. In this case, / is a function of 
two variables u\ and u 2 known through its series expansion in powers of u\ and u 2 - 
The resummation technique used in [5] consists in considering / as a function of U\ and 

z = u 2 /ui\ 

n 

and in resumming with respect to a single variable u\ only. An important hypothesis 
underlying this procedure is that one can safely resum (fTTj) with respect to u\ while 
keeping z fixed, i.e. without resumming with respect to u 2 . Under this hypothesis the 
resummed expression associated with / reads: 

[u(uit; z)} n 



f R ( Ul , z) = V d n (a, a(z), b; z) [°°dt e~* t \ [w{u f ]Z j\ (12) 
V Jo [1 - u(uit; z)\ 



with: 



y/l + a(z)u-l 

ulu; z) = v - = (13) 

y/l + a(z)u + l 

where, as above, the coefficients d n (a,a(z),b,z) in (fl2l) are computed so that the re- 
expansion of the r.h.s. of (fl2"j) in powers of u coincides with that of (fTTj) . 

3. Frustrated magnets. 

The Hamiltonian relevant for frustrated systems is given by: 

n = J d d x{i [{dfaf + (# 2 ) 2 + m\<j>l + ft] + 



(14) 
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where the <f>^ i = 1,2, are A^-component vector fields and u\ and 112 are the coupling 
constants that satisfy U\ > and U2 < 4wj — which corresponds to an Hamiltonian 
bounded from below. For m 2 > the ground state of Hamiltonian ( IT4"1) is given by 
4>\ = 4*2 = while for m 2 < it is given by a configuration where 4>\ an d (j>2 are 
orthogonal with the same norm. The Hamiltonian (114j) thus describes a symmetry 
breaking scheme between a disordered and an ordered phase where the O(N) group of 
rotation is broken down to 0(N — 2) which is relevant for frustrated magnets (for details 
see pQ for instance). 

Let us first recall the FP structure of Hamiltonian (TH|) around d = 4 at leading 
order in e = 4 — d [HI [191 ED]- For iV larger than a critical value N c (d) depending on 
the dimension, the RG equations display, apart the usual Gaussian (it* = u\ = 0) and 
0{2N) (u\ + 0,u* 2 = 0) FPs, two non-trivial (w* ^ and + 0) FPs : one, C+, is 
stable; the other one, C_, is unstable. Above N c (d), the transition is thus predicted to 
be of second order. As A" is lowered starting from values of A^ > N c (d), the two FPs 
C + and C_ get closer and finally collapse together for A" = N c (d). Below N c (d), there 
is no longer a stable FP and the transition is expected to be of first order. The value of 
N c (d) for d = 3 has been computed by several means: in a double expansion in u±, U2 
and e = 4 — d up to five-loops [21 [31 H] , directly in d = 3 in a weak-coupling expansion 
within the massive scheme up to six- loops [6], within a NPRG approach [1] and, finally, 
within the MS scheme without e-expansion [5]. The predictions obtained within the 
perturbative approaches performed at fixed d strongly differ from those obtained using 
other methods, in particular in d < 3, see below. This is the reason which led us to 
reconsider this kind of approach. 

We thus apply the resummation procedure described above without e-expansion 
[2T] to the (3 Ui functions, i = 1,2, obtained at five loops in the MS scheme [5J. More 
precisely, as in [5], we resum (fl Ui (.ui, z ) + £Ui)/ui 2 , 2 = 1,2, instead of (3 Ui (ui, z), which, 
in fact, leads to similar results. For the model ( TT4|) . the region of Borel-summability is 
given by 2u\ — m 2 > (see for instance [5] for details) to which corresponds the value 
a(z) = 1/2. For 2u± — U2 < there exists a singularity on the positive real axis so that 
the series is no longer Borel-summable. However, as noted in [5j, as far as 4ui — U2 > 0, 
the singularity of the Borel transform closest to the origin is still on the negative axis. 
Thus, the asymptotic behavior is still correctly taken into account by the conformal 
mapping and one can, a priori, trust the resummed results. Note finally that b and a 
are typically varied in the ranges [6, 30] and [—0.5, 2]. 

One finds, in agreement with [5], that there exists a curve (parametrized by N c (d) 
or its reciprocal d c (N)) such that for d < d c (N) a stable FP C + governs the critical 
properties of the system. The curves N c (d) obtained within this scheme is shown in FigJT] 
by a dotted line (A^ C FD ). On the same figure we show the results for N c (d) obtained within 
the NPRG approach [I], red solid curve (N^ PRG ) and by (resummed) e 5 -expansion 
[5], green solid curve (Aff). Two points must be noted. First the curves N^(d) and 
N^ PRG (d) show a remarkable agreement being given the very different nature of the 
two corresponding computations. Second, one can see on Fig{T]the strong discrepancy 
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between the two previous approaches and the perturbative FD approach. In particular 
the S-like shape of the curve Nj D (d) obtained within the perturbative FD approach is 
such that the FP C+ D exists for all iV > 2 in d = 3 contrary to the other approaches in 
which a FP C + exists only for N > N c (d = 3) ~ 5 pQ. 

This situation, put together with the fact, already underlined in the Introduction, 
that the FD approach a priori displays spurious FPs, lead us to strongly question the 
validity of the results obtained at FD and, in particular, the existence of genuine FPs 
in d = 3 for iV = 2, 3. 




Figure 1. Curves N c (d) obtained within the MS scheme with e-expansion (N§), 
without e-expansion (iV c FD ) and the NPRG approach (iV c NPRG ). The resummation 
parameters for the MS curve are a = 1/2, b = 10 and a = 1. 

4. Convergence of the loop expansion. 

4-1. The frustrated model. 

We first examine the convergence of the loop expansion of the FD approach by studying 
the sensitivity of the resummed quantities with respect to variations of the resummation 
parameters b and a as well as the order L of the computation. 

We focus on the (real part of the) correction to scaling exponent uj at the FP C + 
which rules its stability : for Re(aj) > the FP is stable and for Re(u;) < it is unstable. 
Within the FD approach, one finds in d = 3 for N < 7 that C + is a focus, that is u is 
complex at this FP (the flow spirals around it). In practice, following [17], we optimize 
u(a, b, L) by first choosing a such that u(L + l)— uj(L) is minimal. Then, one determines 
the parameter b in such a way that u(b) is stationary. We have checked that similar 
convergence properties are obtained for the exponent v |22j. 
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We start by considering the uncontroversial case of a "large" number of components 
of the order parameter, typically greater than 7 in d — 3. Indeed, in this case, all 
perturbative and non-perturbative approaches agree together and find a stable FP 
characterizing a second order behavior pQ. In Figj2j we display u(b) in the case N = 9 
for L = 4 and L = 5 for typical values of a. At a given order, one sees that it is indeed 
possible to find values of b that make u stationary. By performing the same analysis at 
four- and five-loop orders one observes (i) that the dependence of uj with respect to b 
decreases with the order of the expansion as expected (ii) a convergence of the results 
with the order with, however, large error bars typically around 5 — 7%. Note that the 
typical error bars obtained for the Ising model with the same methodology at the same 
orders are much less than 1%. 

CO 

0.8 r 




5 10 15 20 

Figure 2. The critical exponent u> for N = 9 as a function of b at four (curves on the 
left) and five (curves on the right) loops for a = —0.5, 0, 0.5 for the frustrated model. 

We then consider the controversial case of Heisenberg systems (N = 3). In Figj3] 
we display again u(b) for L = 4 and L = 5. There, a new phenomenon occurs: while 
one still finds a stationary value of u(b) at four loops, this is no longer the case at 
five loops. Moreover, considering the values of a and b that minimize the difference 
ui(L = 5) — oj(L = 4) one observes a bad "convergence": the errors bar on the critical 
exponents are now of order 30% and thus far larger than in the case N = 9. An even 
worse behavior is obtained in the XY case. 

From these analyses one gets striking indications of drastically different convergence 
properties for the N = 9 and N = 3 cases. This suggests a qualitative difference between 
the corresponding FPs. To characterize more precisely this difference let us now study 
the cubic model along the same lines. 
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Figure 3. The (real part of the) critical exponent uj for N = 3 as a function of b at 
five (upper curves) and four (lower curves) loops for a — —0.5, 0, 0.5 for the frustrated 
model. Note the change of scale with respect to Fig{2j 



4-2. The cubic model. 

We now consider the ferromagnetic model with cubic anisotropy whose Hamiltonian is: 



with a iV-component vector field 0. The Hamiltonian ( IToT) is used to study the 
critical behavior of numerous magnetic and ferroelectric systems with appropriate order 
parameter symmetry (see e.g. [23]). The /3-functions are known at five- loop order in 
the MS scheme [21] and at six-loop order in the massive scheme [25]. Apart from the 
Gaussian (u* = v* = 0) and an Ising FP (u* = 0,v* ^ 0), there exist two FPs more: 
the O(N) symmetric FP (u* ^ 0,v* = 0) and the mixed one M (u* ^ 0,v* ^ 0). The 
Ising and Gaussian FPs are both unstable for all values of N. The O(N) FP is stable 
and M is unstable with v * < for N < N c and the opposite for iV > iV c . The critical 
value N c has been found to be slightly less than 3: for instance N c ~ 2.89(4) in [25] and 
N c ~ 2.862(5) in [23]. 

Let us now analyze the FP structure of the model (fT5l) within the MS scheme 
without e-expansion by applying the conformal mapping Borel transform ( fl~2l) at d — 3 
(e = 1). The parameter a(z = v/u) entering in (JT2l) is now given by a(z) = 1 + z for 
z > and a(z) = I + z/N for z < while the region of Borel-summability is given by 
the condition u + v > and Nu + v > 0. Within this scheme one surprisingly observes 
that, in addition to the above mentioned usual FPs, there exist, in a whole domain of 
parameters b and a, several other FPs that have no counterparts in the e-expansion. In 
particular, one of them that we call P (which is stable and such that u* > 0, v * < 0) 
exists for any value of iV < 7.5 and lies in the region of Borel-summability u + v > 0. 
The presence of this FP, if taken seriously, would have important physical consequences 




(15) 
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since it would correspond to a second order phase transition with a new universality 
class. However no such transition has ever been reported. On the contrary, a first order 
behaviour for all values of N larger than N c is found within perturbative [25j EE] or 
non-perturbative [27] field theoretical analysis as well as numerical simulations [28] in 
related systems (four-state antiferromagnetic Potts model). Thus, the existence of P 
has to be considered as an artefact of the FD analysis. Note finally, and interestingly, 
that P is found to be a focus FP, a striking similarity with the frustrated case. 




At this stage, it is very instructive to perform the same convergence analysis as the 
one performed for C + in the frustrated case. The Ising, O(N) and mixed M FPs display 
good convergence properties and we focus in the following on the supernumerary FP 
P. In FigH] we plot u(b) at this FP for L = 4, 5 and for three different values of a. 
Interestingly, when comparing Figj3]and FigJUwe find similar behavior between the cubic 
case and the frustrated case for N = 3. Indeed, on FigJUone finds stationary values of 
u(b) at four loops but, at five loops, this is only the case for a = 1.7. It is also interesting 
to consider the values of a and b that minimizes the difference u(L = 5) — u(L = 4). 
From there, one finds error bars for u of order 40%, i.e. of the same order of magnitude 
as the one found in the N = 3 frustrated case. Being given the fact that the FP P is 
clearly an artefact of the FD analysis, our study suggests that the lack of convergence of 
ui{b) characterizes the behavior at a spurious FP. Thus, coming back on the frustrated 
case one is naturally led to the conclusion that the FPs C+ D found in d — 3 and iV = 2, 3 
should also be interpreted as spurious FPs, as artefacts of the FD analysis. In order to 
confirm this statement we now examine another characteristic feature of the set of CY° 
FPs considered as functions of d and N. 
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5. The singularity S. 

We now display a specific feature of the RG flow of frustrated magnets analyzed within 
the FD approach that provides another strong indication of the problematic character 
of this approach. Let us consider the coordinates (u*,u 2 ) of the FP C+ D as functions 
of d and N: u* = u*(d, N), i = 1,2. These functions are the roots of the /3-functions of 
the couplings U\ and u 2 obtained in the FD approach and resummed according to the 
scheme sketched above, Eq.(fT2l): 

f3 Ul (ul,u* 2 ) = f3 U2 (ul,u* 2 ) = . (16) 

The resummed /3-functions are smooth functions of d and N showing no particular 
feature for 2 < d < 4 and 2 < iV < oo. However, as we now show, the functions 
it* = u*(d, N), i = 1, 2 exhibit a non-trivial behavior as d and N are varied continuously 
around the point labelled S in FigUJ 

Let us first give, in FigJHl the precise definition of the curve N c (d) in the FD 
approach. This curve is made of two parts: the first one, labelled (I), corresponds to 
the part of the curve above S while the second one, labelled (II), corresponds to the 
part below S. (I) is defined as the line in the (d, N) plane above which there exist two 
non-trivial locus FP (that is having real exponents uj), one stable C+ D and one unstable 
CY°, and below which there exists none. We recall the mechanism of disappearance 
of these FPs: when at fixed dimension d, N is decreased from large values down to 
iV = N c (d) the two FPs and CYP get closer and closer from each other and finally 
collapse right on (I). For N below N c (d), the coordinates u\ and u 2 become complex 
and there is no longer any non-trivial FP with real coordinates. Thus, the part (I) of 
the curve N c (d) corresponds to the region where the speed of the RG flow between CY° 
and CYP vanishes. The same behavior is observed in other (e-expansion and NPRG) 
approaches and the numerical values of N c (d) in the corresponding part of the curve are 
very close when calculated by different approches (c.f. FigH]). 

As usual, the speed of the RG flow around a FP can be obtained by linearizing the 
flow equations at this point. It is thus governed by the two eigenvalues of the "stability" 
matrix 

d(3 Ui (ui,u 2 ) 



;i7) 

**3 u^,v,2 

that provide the speeds of the RG flow in its two eigendirections at the FP considered. 
These eigenvalues define the two critical exponents u)\ and uo 2 governing the corrections 
to scaling. The equation of the part (I) of N c (d) is thus given by: 

uj\ = oj\ (u*(d, N),u 2 (d, N)\ = (18) 

where U\ is the eigenvalue of M corresponding to the eigendirection of the flow joining 
C™ to CYP ■ Now, when moving on the line N c (d) towards S one finds that u 2 decreases 
and eventually vanishes. One thus defines the point S by: 

lu 2 (uI(S),u* 2 (S)) = 0. (19) 




Figure 5. Curve N c (d) within the MS scheme without e-expansion. The part (I) of the 
curve iVj corresponds to a boundary between a region where, at a given dimension d 
there exists a stable locus FP at large N and no FP at small N. The part (II) of the 
curve iVj D corresponds to a boundary between a region where at fixed N there exists 
an attractive focus FP for d < d(N) and a repulsive focus for d > d(N). Finally the 
line (F) is defined as follow: above (F), C+ D is a locus FP (u>\ and 0J2 are both real) 
and below (F), it is a focus (u>i and L02 = u>* are complex). 

Thus, right at S, uj\ = uj 2 = and, with the choice of parameters a = 1/2, b = 10, a = 1 
one finds that at 3, d = 3.24, N = 7 and K(S), u* 2 (S)) = (0.3,0.7). Thus, S is just 
outside — but not far of — the region of Borel-summability. Note also that, since its 
coordinates verify 4m* — u 2 > 0, S is still in the region where the resummed results are 
trustable. 

In fact, the point S is also a special point in the sense that it is the endpoint of 
another line, labelled (F) on Fig 15] which is such that below (F), the two exponents 
Wj's acquire an imaginary part and are complex conjugates: U\ = This means that, 
below (F), the FPs become focuses, either attractive or repulsive. This is in particular 
the case of the FP CY° found in [3] in d — 3. From now on we concentrate on the 
FP C™ since the coordinates of CYP rapidly grow and go outside the region of Borel 
summability. Note that the part of (F) shown on Figj5] lies inside the region of Borel 
summability and thus, within the FD approach, is supposed to be well under control. 

We now define the part (II) of the curve N c (d) as the line separating the region 
where C+ D is an attractive focus FP (for d < d c (N)) and the region where it is a repulsive 
focus FP (for d > d c (N)) 0. Thus, by definition, (II) is the line on which the real part 

\ The change of stability of C_j_ D occurs as follows. For d < d c (N) and N < 7.5 and sufficiently close to 
(II), C™ is a focus and is attractive inside a basin of attraction which is a limit cycle for the RG flow. 
This limit cycle is repulsive : inside it the RG flow converges to CY°, outside the RG flow diverges. 
When, at fixed N, the dimension d is increased this limit cycle shrinks and becomes a point on (II). 
When d is further increased, the limit cycle re-appears but now it is attractive and C+ D becomes a 
repulsive focus FP. 
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of the Wj's vanishes: 

Re(wi) = Re{u 2 ) = (20) 

but it is no longer a line where Ci° collapses with another FP and disappears. It is 
remarkable that within the FD approach, there exists a line (the part (II) of the curve 
N c (d)) that has no counterpart in other approches. 

Contrary to what occurs above, the coordinates u\ and u\ of C+ D on the line 
(II) are, for a large part of (II), far outside the region of Borel summability. It is 
thus not possible to compute accurately the location of (II). We, however, emphasize 
that only the existence of (II) is necessary for the validity of our arguments, not its 
precise determination. As for the existence of this part (II) of the curve N c (d), it is an 
unavoidable consequence of the existence of S which is supposed to be under control 
within the perturbative FD approach. Thus, either S really exists and part (II) of N c (d) 
also exists — with a shape that could be somewhat different of the one drawn on Figj5] 
- or it does not exist and neither does S. In this second case, this would mean that the 
whole resummation scheme is questionable at least for sufficiently low d and iV (typically 
d < 3.2 and N < 7). We argue in the following that this is very probably what occurs. 

Let us now show on Figj6] that there exists a very non-trivial property of the RG 
flow which is a consequence of the existence of S. The idea is to follow continuously the 
coordinates of the FP C+ D along a path encircling S, path A in Figj6]for instance. We 
start, for instance at (d = 3, N — 5) go to (d = 3.4, N = 5) then to (d = 3.4, N = 9) 
then to {d = 3,N = 9) and finally go back to (d = 3,N = 5). The surprising fact is 
that by making a trip along such a closed path, the coordinates {u\,u^} of C? D do not 
go back to their original value. This is specific to S since along any closed path that 
does not encircle this point — path B in Figj6]for instance — the coordinates {u\,u^) 
of C+ D always go back to their original values. Let us emphasize here that such a path 
can well cross part (I) of the curve N c (d) as path B does in Figj6j In this case, C+ D has 
complex coordinates on the part of the path which is below (I), but as the path crosses 
again (I) the coordinates become real again and go back finally to their original values. 
All this makes S a topological singularity of the functions u\(d, N) and u^(d, N). 

Let us remark at this stage that, before having drawn any conclusion from 
the existence of a topological singularity, one faces a very unusual property of the 
mapping from the plane (N, d) to the FP coupling constants space (u*, u*£) which makes 
multivalued the functions u*(N, d) and u^(N, d). While one cannot a priori discard the 
possibility of such a behavior of the functions ul(N,d) and u^N, d), it is tempting to 
attribute it to the coexistence at a given (d, N) of FPs that are identical to those found 
within the e = 4 — d expansion and FPs that are artefacts of the FD approach. Let us 
now draw the full consequences of the striking behavior described above. 

6. Fixed point at the upper critical dimension. 

We now present our last argument in favor of the spurious character of the FPs obtained 
within the FD approach for small d and N. It is based, to a large extent, on the existence 
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3 3.2 3.4 3.6 3.8 4 

Figure 6. Two different paths in the (N, d) plane. Path A encloses S and is such 
that the coordinates of the FP CT D do not go back to their original values after a trip 
along this path at variance with path B. 



of the singularity S in the (u*,v%) plane that leads to a striking property of the field 
theory describing frustrated magnets. 

Let us first recall some basic features about the description of the long distance 
physics of a lattice system by a continuous field theory. The most important ingredient in 
this construction is the choice of a low-energy effective Hamiltonian in terms of the order 
parameter 0. This choice implies the selection of a finite number of terms — the 2 and 
4 terms for second order phase transitions — among the infinite number of operators 
obtained in the Hubbard-Stratonovitch derivation of the microscopic hamiltonian. This 
selection of the most relevant terms fully relies on the existence of an upper critical 
dimension where the theory is perturbatively infrared free, i. e. controlled by the Gaussian 
FP. Indeed, it is only under this condition that power-counting makes sense since it is 
based on the engineering dimension of the field that neglects fluctuations. Perturbation 
theory can then be used since, by definition, it consists in an expansion around the 
Gaussian theory. Also, the perturbative results obtained in this way are reliable as long 
as the (infrared attractive) FP found this way (i) is connected to the Gaussian by the 
RG flow, (ii) lies in the Borel-summability region (iii) is not too far from the Gaussian 
FP so that calculations performed at L loops lead to converged results. Let us indicate 
that because of point (i) above, if the theory is trivial in d = 4, it is very probable that 
any non-trivial FP identified in d — 3, once followed continuously from d — 3 to d — 4, 
becomes Gaussian in this dimension. This is in particular what underlies the validity of 
the e = 4 — d expansion. 

Let us now briefly discuss the question of triviality of scalar theories in d — 4 
and the ensuing consequences for their perturbative analysis in any dimension below 4. 
We first emphasize that there is no rigorous proof of the triviality of scalar theories in 
d = 4. However (i) there is a large body of evidences of triviality at least for the O(N) 
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models and in particular for the Ising model (ii) it is very likely that even if the scalar 
models were non-trivial in d — 4, perturbation theories would not be able to reach the 
corresponding FP [29J. Thus, the most natural hypothesis is that the O(N) x 0(2) 
theory is also trivial in d = 4. An ever weaker hypothesis, that we make and use in 
the following, is that this is the case when it is analyzed perturbatively. All previous 
considerations and assumptions lead us to conclude that, very probably, any physical 
FP found perturbatively in any dimension must be Gaussian when followed continuously 
in d = 4. Let us notice that a non-perturbative approach, performed by some of the 
present authors [1] on the O(N) x 0(2) theory did not lead to any non-trivial FP in 
d — 4 which sustains our triviality hypotheses in d — 4. 

Therefore, according to our hypotheses, a practical way to check whether a FP 
found at a given dimension, d = 3 for instance, is a genuine FP or is just an artefact of 
perturbation theory is to follow it by continuity up to d = 4 [3j [30] . If the FP survives 
as a non-Gaussian FP at this dimension we consider it as spurious. It is important 
to realize that our criterion does not exclude FP in d = 3 whose coordinates become 
complex between d = 3 to d = 4 as far as they vanish in d — 4. This is in particular the 
case for all FPs associated to the paths with fixed N in the (d, N) plane that cross (I) 
above S. 

Let us apply this criterion to analyze the FPs P and CT D that appear in the FD 
analysis of the cubic and frustrated models. We present our results in FigJT] where we 
have displayed the coordinates u* and u\ associated with these FPs as functions of d at 
fixed N. Manifestly, they both survive everywhere above d = 3 and are not Gaussian 
in d = 4. In the cubic case, this is true for all values of N for which P exists. In the 
frustrated case, this is true for N typically below 7. According to our criterion, P is, as 
expected, always found to be spurious while C+ D is spurious only for N < 7. We thus 
conclude, in the frustrated case, that the FPs found in d = 3 for N = 2, 3 are spurious. 



* * 
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Figure 7. The u* coordinate of the FP P of the cubic model (N = 2, upper curve) 
and the u\ coordinate of the FP C™ of the frustrated model (N = 3, lower curve) as 
functions of d. 
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Let us emphasize that when d approaches 4 the coordinates of the FPs P and CT D 
for N < 7 become large and do no longer belong to the region of Borel-summability 
so that they cannot be determined accurately as it is the case of part (II) of the curve 
N c (d). One could thus naively conclude that it is not possible to safely decide whether 
these FPs are or are not Gaussian in d — 4. This is actually not the case for at least 
two reasons. First, if these FPs were Gaussian in d — 4, their coordinates just below 
this dimension would be extremely small and their existence as Gaussian FPs in d — 4 
could be safely established within perturbation theory even without any resummation 
procedure. Since for d just below 4 no such FPs close to the Gaussian are found in 
perturbation theory, their non-Gaussian character is doubtless in d — 4. Second, in the 
frustrated case, the non-Gaussian character of the C™ FP is a clear consequence of 
the existence of (II) which, itself, relies on the existence of the singularity S. Indeed, 
following a FP C+ D along a path starting in d = 3, going to d = 4 and crossing Nj D (d) 
above S, the coordinates of become complex in d = d c (N) and both the real and 
the imaginary parts of u\ and u\ go to zero for d = 4 where it is thus a Gaussian FP. If, 
on the contrary, the path crosses Nj D (d) below S, CT° changes from a stable focus to 
an unstable one at d = d c (N) and does not go to the Gaussian in d — 4. The singularity 
S is therefore responsible for the non-Gaussian character of CT D in d = 4. Thus, even 
if the coordinates of C+ D for N < 7 become large for d close to 4, its non-Gaussian 
character in this dimension makes no doubt. Note finally that it is tempting to follow 
this FP CT° above d = 4 where there exist rigorous proofs of the triviality of the scalar 
4 field theory [3 lj - Indeed, as it is suggested by FigfTl the FP C+ D apparently survives 
at finite distance above d=4. However this last fact must be taken with great caution 
since C+ D is now deeply out of the region of Borel-summability and one has no longer 
any control of where the FP really lies. 

7. Conclusion. 

It appears from our study that the FPs identified in the FD approach are very 
likely spurious. The transition in frustrated magnets should thus be of — possibly weak 
- first order in agreement with NPRG and e-expansion approaches. It remains to 
explain the failure in the resummation procedure used in the FD approach, Eq. (fT2]) . As 
already emphasized this procedure relies on the hypothesis that resumming with respect 
to Mi, keeping a polynomial structure in w 2 , is sufficient. Alternatively, a resummation 
of the series with respect to the two coupling constants could be required to obtain 
reliable results (see for instance [32] for the randomly diluted Ising model). Postponing 
these considerations for a future publication [22J we assume that the use of Eq. ( Tl2l) as 
such is justified. Then, a possible origin of the failure in the resummation procedure 
could be that the series considered are not known at large enough order to reach the 
asymptotic behavior. In this case there would be no reason to fix the parameter a at 
its asymptotic value a = 1/2 and one should have to vary it as b and a to optimize the 
results [TTj . We display in FigJHlthe curves Nf B (d) for different values of a. The part 
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corresponding to large values of N c , typically for iV c > 7, is almost insensitive to the 
variations of a whereas this is clearly not the case for smaller values of N c . In particular, 
for sufficiently large values of a, typically a > 1.5, the S-like part is pushed below d = 3 
so that N^(d), Nj B (d) and N^ PRG (d) are then compatible everywhere for 3 < d < 4. 
Let us also notice in this respect that, for a = 1.3, the shape of the curve N™(d) is 
even compatible with the results obtained in the massive scheme in d = 3 in which one 
finds FPs for all values of N but in the range 5.7(3) < N < 6.4(4) [33J. This suggests 
that the two FD methods — massive scheme in d = 3 and MS without e-expansion - 
are in fact compatible but suffer from the same problems of convergence. Thus, under 
our hypothesis, all qualitative differences between the different approaches disappear, 
so that the problem would boil down to a question of order of computation. 



Nc 




2 . 6 



d 



Figure 8. Three curves N^ D (d) for different values of the parameter a (from right to 
left a = 0.5, 1.3, 1.5) and the curve N^(d). The parameters b and a are b = 10 and 
a = 1. The parts of the curves below the black dots correspond to a regime of Borel 
non-summability. 



Finally note that our present considerations surely pertains to the case of frustrated 
magnets in d = 2 [34] . Indeed we have checked that the FP found in d = 2 is continuously 
related to the FP C+ D in d = 3 which makes its existence doubtful. Our conclusions 
could also apply in other situations where FPs that have no counterpart in the e- 
expansion are found, as it is the case, for instance, in QCD at finite temperature [35] . 
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